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Infrared singularities of scattering amplitudes in perturbative QCD 
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An exact formula is derived for the infrared singularities of dimensionally regularized scattering am- 
plitudes in massless QCD with an arbitrary number of legs, valid at any number of loops. It is based 
on the conjecture that the anomalous-dimension matrix of n-jet operators in soft-collinear effective 
theory contains only a single non-trivial color structure, whose coefficient is the cusp anomalous 
dimension of Wilson loops with light-like segments. Its color-diagonal part is characterized by two 
anomalous dimensions, which are extracted to three-loop order from known perturbative results for 
the quark and gluon form factors. This allows us to predict the three- loop coefficients of all 1/e 
poles for an arbitrary n-parton scattering amplitudes, generalizing existing two-loop results. 



Perturbative expressions for quark-gluon scattering 
amplitudes in massless QCD contain infrared (IR) singu- 
larities, which originate from configurations where loop 
momenta become soft or collinear. In predictions for 
physical observables these cancel against corresponding 
singularities from real gluon emission. Nevertheless, the 
IR singularities are interesting in their own right. Af- 
ter the cancellation of IR singularities logarithmic terms 
remain, which depend on the phase-space cuts imposed. 
For high-energy scattering processes with low-mass jets 
in the final states these Sudakov logarithms dominate the 
cross section. An understanding of the structure of IR 
singularities can be used to predict and resum these log- 
arithmically enhanced contributions to all orders. It also 
serves as a consistency check on loop calculations. 

Despite the fact that their origin is well understood, an 
all-order description of the singularity structure of QCD 
amplitudes is currently still lacking. An important step 
toward this goal was made by Catani [J], who correctly 
predicted the singularities of two-loop amplitudes apart 
from the 1/e pole term. His formula states that the prod- 
uct 
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where a s = a s (p), and A4 n (e, {p})) with {p} = 
{px, . . . ,p n } denotes an ultraviolet (UV) renormalized, 
on-shell n-parton scattering amplitude with IR singu- 
larities regularized in d = 4 — 2e dimensions, is free 
of IR poles through 0{a s ). The subtraction operators 
J(«)( e ) = J(«)(e ; {p}, ,j) are defined as 
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where Sij = 2o~ij pi -pj +i0, and the sign factor <7y = +1 
if the momenta pi and pj are both incoming or outgo- 
ing, and aij = —1 otherwise. The quantities K and 
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([2]) are related to anomalous-dimension coefficients 
below as K = r™ sp /(2r™ sp ) and g, = - 7 §/2, 
while /3q = Ca — § ?i?n/ is the first coefficient of the 
QCD /3-function. The sums in the expressions above are 
over all external partons. We use the color-space for- 
malism of , in which n-particle amplitudes are treated 
as n-dimensional vectors in color space. T\ is the color 
generator associated with the «-th parton and acts as 
a matrix on its color index. Specifically, one assigns 
(Ti) a p — ta8 f° r a final-state quark or initial-state anti- 
quark, {T°~) a p = —t a p a for a final-state anti-quark or 
initial-state quark, and (T")i, c = —if a bc for a gluon. The 
product Tj • Tj = T° is summed over a. Genera- 
tors associated with different particles trivially commute, 
Ti ■ Tj = Tj ■ T t for i^j, while T? = C l is given in terms 
of the quadratic Casimir operator of the corresponding 
color representation, i.e., C q = Cq = Cf and C g = Ca- 
Owing to color conservation, the scattering amplitudes 
fulfill the relation J2i T i |-M„(e, {p})) = 0. Because the 
operator T, ■ Tj is color-conserving, this condition is also 
fulfilled after acting with or on the amplitudes. 

(2) 

The scheme-dependent quantity s , which only 
contains 1/e divergences, was not specified in [l[ except 
for the simplest case of the on-shell quark form factor. In 
subsequent two-loop calculations of the three-parton qqg 
amplitude [|[ and of 2 — ► 2 scattering amplitudes 0, 0] it 
was observed that its color-diagonal part has a universal 
structure. However, the calculations of four-parton am- 
plitudes also revealed the presence of a non-trivial color 
structure Q. A conjecture for the form of this term for 
a general n-parton amplitude was made in Q . 

An interesting alternative approach to the problem of 
IR singularities of on-shell amplitudes was developed in 
[7[ , where the authors exploited the factorization proper- 
ties of scattering amplitudes [1, Q along with IR evolu- 
tion equations familiar from the analysis of the Sudakov 
form factor [lfj. They recovered Catani's result (J2J) at 
two-loop order and related the coefficient of the unspeci- 
fied 1/e pole term to a soft anomalous-dimension matrix, 
which was unknown at the time. They also explained how 
their method could be extended beyond two-loop order. 
The two-loop soft anomalous-dimension matrix was later 
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calculated in where its color structure was found to 
be proportional to that obtained at one- loop order. 

In this Letter we propose an all-order generalization of 
Catani's result @ valid for an arbitrary on-shell n-parton 
scattering amplitude. We find that in a minimal subtrac- 
tion scheme the color structure of the IR pole terms is 
simpler than anticipated based on Catani's work [l[. In 
fact, to all-loop order the 1/e pole term contains only 
the structures 1 and Ti • Tj. Our analysis is based on 
effective field theory and shares many similarities with 
that of (7j. However, in our case the hard, jet, and soft 
functions are defined in terms of matrix elements of dif- 
ferent types of fields in the effective theory and are in 
one-to-one correspondence with different physical scales. 
The corresponding definitions in Q are less intuitive. 

Our key observation is that the IR singularities of on- 
shell amplitudes in massless QCD are in one-to-one corre- 
spondence to the UV poles of operator matrix elements in 
soft-collinear effective theory (SCET) [11 [3. They can 
be subtracted by means of a multiplicative renormaliza- 
tion factor Z (a matrix in color space), whose structure 
is constrained by the renormalization group (RG). SCET 
is the appropriate effective theory to analyze scattering 
processes at large momentum transfer, which involve jets 
(or individual hadrons) with small invariant masses. It 
separates hard contributions associated with the large 
momentum transfer from low-energy contributions as- 
sociated with the small invariant masses of the initial- 
and final-state particles. For a general n-jet observable, 
the effective theory involves a set of collinear fields for 
each direction of large energy flow, which describe the 
QCD dynamics inside the individual jets. It also con- 
tains soft quark and gluon fields, which mediate low- 
energy interactions among the jets. Hard interactions 
are integrated out and absorbed into the Wilson coeffi- 
cients of operators built from soft and collinear fields. A 
generic n-jet process is mediated by an effective Hamilto- 
nian 7i. n = X^i ^n,i{^) 0£™(u), where the sum runs over 
a basis of SCET operators built from n distinct types of 
collinear fields. The bare matrix elements of these oper- 
ators are UV divergent and are renormalized in the MS 
scheme. Their divergences are absorbed into a renormal- 
ization factor via 0, r ™(") = ^ Zytp, e) Oj$°(e). For 
physical quantities, the scale dependence of the Wilson 
coefficients C„.i(u) cancels against that of the matrix el- 
ements of the renormalized operators. 

In a physical process with initial- and final-state 
hadrons, the soft and collinear scales are set by nonper- 
turbative dynamics or experimental cuts. Let us now 
consider (slightly) off-shell n-parton amputated Green's 
functions G n ({p}). In this case the jet-scale Aj is set 
by the off-shellness p\ of the fields, and the soft scale is 
A s ~ Aj/Q, where Q is a typical hard momentum trans- 
fer. The Green's functions are obtained by taking matrix 
elements of the above effective Hamiltonian, which can 
be written as 

G„(M) = lira C «» z a(M, e ) (<%T(e)) . (3) 
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where we suppress the dependence of the quantities on 
the right-hand side on the parton momenta. To ob- 
tain on-shell n-parton scattering amplitudes from these 
Green's functions one takes the limit p\ — > 0. This intro- 
duces IR divergences, which can be regulated by evalu- 
ating the effective-theory matrix elements in d = 4 — 2e 
dimensions. Doing so renders the matrix elements of the 
operators 0^ aro trivial: in the limit pf — > both the soft 
and the jet scales tend to zero, and all loop diagrams 
in the effective theory become scaleless and vanish. The 
bare matrix elements arc thus reduced to trivial Dirac 
and color structures. Since the IR divergences are inde- 
pendent of the spin structure, we will not make the Dirac 
structures explicit but simply absorb them into the Wil- 
son coefficients. The on-shell Green's functions are then 
directly proportional to the Wilson coefficients of n-jet 
SCET operators in the MS scheme. In the color-space 
basis notation of @, the effective Hamiltonian reads 
7~L n = (0™ n |C„), and we have 

|C„(M,M)> = lim n Z- 1 (e,{p}^)\G n (e,{p})). (4) 

This notation is convenient but unconventional, in that 
our Wilson coefficients and operators are not separately 
color singlets and Lorentz scalars. The scattering ampli- 
tudes |A4„(e, {p})) are obtained by contracting the am- 
putated on-shell Green's functions with the spinors and 
polarization vectors associated with the external parti- 
cles. Their singularities are thus governed by the same 
Z matrix. 

The logarithm of the renormalization factor Z in 
is related via T = —dhxZ/dhifi to the anomalous- 
dimension matrix T governing the RG evolution of the 
n-jet SCET operators 0^ n . The same quantity controls 
the evolution of the Wilson coefficients, and hence of the 
minimally subtracted on-shell scattering amplitudes, via 
the evolution equation 

-^—\C n {{p},n)) = T\C n {{p}, l i)). (5) 
din u 

We will now present a conjecture for the exact form of 
the anomalous-dimension matrix. In general, T = T c+S 
is determined by the anomalous-dimension contributions 
of collinear and soft modes in the SCET matrix ele- 
ments. An important feature of SCET is that the in- 
teractions of collinear fields with soft gluons can be re- 
moved by field redefinitions and absorbed into soft Wil- 
son lines [T3 |. Interactions with soft quarks are power 
suppressed and can be ignored. Moreover, the different 
collinear sectors in SCET do not interact with each other. 
This allow us to decompose T = F s + X)i7cj where the 
one-particle collinear contributions are diagonal in color 
space. Hence, contributions to the anomalous dimension 
involving correlations between several partons only reside 
in the soft contribution T s . After the decoupling trans- 
formation the soft matrix element is a vacuum expecta- 
tion value (0\Si . . . S n \0) of n light-like Wilson lines, one 
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for each external parton (see also 0]). Here 



Si = P exp 



ig / din, ■ A a (trii)T? 
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where rii is a light-like reference vector aligned with the 
«-th particle's momentum. When the color indices are 
contracted into color-singlet combinations, the soft ma- 
trix element consists of one or more Wilson loops (closed 
at infinity), whose lines cross and or have cusps at the 
origin. The renormalization properties of Wilson loops 
have been studied extensively in the literature (see e.g. 
0, HE EE E3)- Based on these results and on other 
considerations to be explained in detail elsewhere, we 
propose that the exact expression for the anomalous- 
dimension matrix of an n-jet operator in SCET in the 
color-space formalism is 



in d = 4 — 2e dimensions, we find 
da 
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where we have defined 
d 
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When integrating over a in JSJ, the scale a in the argu- 
ment of the logarithm in T must be kept fixed. Note that 
when acting on color-singlet states, the unweighted sum 
over color generators satisfies the relation 



(10) 



r = T * T i r cus P (a s ) In 3— + ^(a.) • (?) 

The sums run over the n external partons. Here and 
below the notation (ii, ...,ik) refers to ordered tuples of 
distinct parton indices. Our result contains three univer- 
sal anomalous-dimension functions. The quantity r cusp is 
proportional to the cusp anomalous dimension of Wilson 
loops with light-like segments [TtJ , while 7 9 = 7* and j 9 
are anomalous dimensions specific for (anti-)quark and 
gluon fields. They are defined by relation ([7]). 

Our conjecture for the anomalous-dimension matrix 
is the simplest expression consistent with existing cal- 
culations. It implies that the cusp anomalous dimen- 
sion characterizes the renormalization of Wilson lines 
even in the general case, where several lines meet at a 
single space-time point. The structure of the logarith- 
mic terms in ([7J agrees with an explicit two-loop cal- 
culation in 11]. For the quark and gluon form factors, 
the divergent terms are known to three-loop accuracy 
fl8l ]. When applied to this case, our general relation 
([7]) implies that the cusp anomalous dimensions in the 
fundamental and adjoint representations of SU (N c ) are 



related by rf usp (a s )/C F 



T^ sp (a s )/C A 



This relation is indeed fulfilled to three-loop order [19j . 
The application to the form factors also determines the 
anomalous dimensions 7 9 = 7^ and 7 s to three-loop accu- 
racy. In the effective theory the form factors are mapped 
onto two-jet operators containing two collinear quark 
or gluon fields. The corresponding three-loop anoma- 
lous dimensions were given in [201 ] . In the notation of 
these papers, we have 2-f q (a s ) — j v (a s ) and 2j 9 (a s ) = 
7*(a s ) + j s (a s ) + (3(a s )/a s , where (3{a s ) = da s /d\nu. 
These results are valid in the standard dimensional reg- 
ularization scheme with d-dimensional helicities. 

The relation T = — dlnZ/dlnfj, may be integrated to 
obtain a closed expression for the logarithm of Z. Using 
the relation /3(a s ,e) = (3(a s ) — 2ea s for the /3-function 



which follows from color conservation. Since both the 
scattering amplitudes and the anomalous dimension T 
are color-conserving, this relation can always be used 
in our case. It is understood that the result © must 
be expanded in powers of a s with e treated as a fixed 
O(a s ) quantity. Writing the perturbative series of the 
anomalous dimension and /^-function in the form T(a s ) — 
E„T„(^)"+ 1 and 0(a s ) = -2a s £„ (3 n {^) n+ \ we 
find up to three-loop order the result 
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which only contains the two color structures present in 
(|7]). Exponentiation yields an explicit expression for Z, 
which simplifies since the different expansion coefficients 
r„ commute. The relevant one-loop coefficients of the 
three anomalous-dimension functions are Tq USP = 4, 7q = 
— 3CV, and 7$ = — (3q. The highest pole in the C?(a") 
term of \nZ is l/e Tl+1 , instead of l/e 2 ™ for the Z factor 
itself. The exponentiation of the higher pole terms was 
noted previously in Q. 

The one- and two-loop coefficients of the matrix Z 
are closely related to Catani's subtraction operators 1^ 
and /( 2 ) given in The conditions linking his objects 
to ours are that the differences 21^' — Z\ and 4J( 2 ) + 
remain finite for e — » 0. Here Z n denotes 
the coefficient of (a s /47r)™ in the Z factor. The first 

relation is indeed satisfied. The second one can be used 

(2) 

to derive an explicit expression for the quantity H\ s 
in @ entering the two-loop coefficient of the l/e pole, 
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which was not obtained in [l[ . We find 

H nl = We £ (ft ~ ff» ^ + is h Cij (12) 

+ ifabc T?T*T£hi^l- ln^ln^i, 
4e • -s jk -s kl s i:j 

which apart from the last term is diagonal in color space 
and universal in the sense that it is a sum over contri- 
butions from individual partons. Note that only the first 
term in this result is of a form suggested by (fTTj) . The 
remaining terms in the first line arise because two-loop 
corrections involving the cusp anomalous dimension or 
the /3-function are not implemented in an optimal way 
in @. More importantly, the last term in (fHi|) . which 
is present for four or more partons, arises only because 
the subtraction operators 1^ in [l[ are not defined in a 
minimal scheme but also include 0(e k ) terms with k > 0. 

Our expressions (fTTj) and (fT^]) reproduce all known re- 
sults for the two-loop l/e k poles of on-shcll scattering 
amplitudes in massless QCD. In addition to the on-shell 
quark and gluon form factors, these include e + e~ — > qqg 
3[ as well as all four-point functions of quarks and gluons 
J, [f|. It further confirms the ansatz made for higher- 
point functions in |6j. At the three- loop level, only the 
IR divergences of the quark and gluon form factors are 
known for the QCD case [l8| . For J\f — 4 supersymmet- 
ric Yang-Mills theory in the planar limit, on the other 



hand, the four-point functions are known up to four-loop 
order [2l| . The divergent part of these amplitudes factors 
into a product of square roots of form factors of neigh- 
bouring legs, which is consistent with the structure of 
our anomalous dimension, given that at leading order in 
1 /N c one has Ti • 7} — > — 4^ for neighbouring legs and 
zero otherwise. 

Let us briefly mention some applications and possi- 
ble generalizations of our results. The most important 
one is the resummation of Sudakov logarithms in n-jct 
processes. The evolution equation ([5]) is simple enough 
to admit exact solutions in closed form, and the known 
three-loop anomalous dimensions allow for resummations 
at next-to-next-to-next-to-leading logarithmic accuracy. 
Our formalism can be generalized to processes involving 
massive partons by combining it with methods developed 
in (221 . [23( 1 . Since the IR poles must cancel between vir- 
tual and real corrections in a given physical process, our 
results for the virtual corrections could help to improve 
the analysis of real-emission corrections as well. More- 
over, the great simplicity of our result ([7]) appears to 
hint at a universal origin of IR singularities that is dis- 
connected from the genuine dynamics of the scattering 
amplitude itself. 
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